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ABSTRACT 

We show that all X-ray decay curves of GRBs measured by Swift can be fitted 
using one or two components both of which have exactly the same functional form 
comprised of an early falling exponential phase followed by a power law decay. 
The 1st component contains the prompt 7-ray emission and the initial X-ray 
decay. The 2nd component appears later, has a much longer duration and is 
present for fh 80% of GRBs. It most likely arises from the external shock which 
eventually develops into the X-ray afterglow. In the remaining « 20% of GRBs 
the initial X-ray decay of the 1st component fades more slowly than the 2nd and 
dominates at late times to form an afterglow but it is not clear what the origin 
of this emission is. 

The temporal decay parameters and 7/X-ray spectral indices derived for 107 
GRBs are compared to the expectations of the standard fireball model including 
a search for possible "jet breaks". For ~ 50% of GRBs the observed afterglow is 
in accord with the model but for the rest the temporal and spectral indices do 
not conform to the expected closure relations and are suggestive of continued, 
late, energy injection. We identify a few possible jet breaks but there are many 
examples where such breaks are predicted but are absent. 

The time, T a , at which the exponential phase of the 2nd component changes 
to a final powerlaw decay afterglow is correlated with the peak of the 7-ray 
spectrum, E pea k. This is analogous to the Ghirlanda relation, indicating that 
this time is in some way related to optically observed break times measured 
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for pie-Swift bursts. Many optical breaks have previously been identified as jet 
breaks but the differences seen between X-ray and optical afterglows suggest that 
this is not the explanation. 

Subject headings: Gamma Rays: bursts — radiation mechanisms: non-thermal 
— ISM: jets and outflows 

1. Introduction 

The standard fireball shock model of GRBs (Meszaros 2002 and references therein) 
predicts that a broadband continuum afterglow spectrum is expected to arise from an ex- 
ternal shock when the relativistically expanding fireball is decelerated by the surrounding 
low density medium. As relativistic electrons, accelerated in the shock to form a power law 
energy spectrum, spiral in the co-moving magnetic field we should see a characteristic fading 
synchroton radiation spectrum stretching from radio frequencies through the IR, optical and 
UV bands in to an X-ray and gamma ray high energy tail. The detailed form of the expected 
afterglow spectrum and its evolution are described by Sari, Piran and Narayan (1998) and 
Wijers and Galama (1999). 

X-ray afterglows of GRBs were first detected by the Beppo-SAX satellite (1996-2002) 
and the detection of GRB970228 (Costa et al. 1997) and other X-ray afterglows provided 
positions of sufficient accuracy to enable follow-up ground-based optical observations. Faint 
optical afterglows were discovered and it was soon established that GRBs occurred at co- 
mological distances. The first redshift, z=0.835, was measured for GRB970508 (Metzger et 
al. 1997). A connection between GRBs and supernovae was revealed by observations of 
GRB980425/SN1998bw (Galama et al. 1998, Kulkarni et al. 1998) although the supernovae 
associated with GRBs showed very high expansion velocities (tens of thousands of kilometers 
per second) and were given a new classification of hypernovae. The XMM-Newton obser- 
vatory also detected X-ray afterglows. In particular GRB030329 confirmed the hypernova 
connection (Tiengo et al. 2003, Stanek et al. 2003, Hjorth et al, 2003) and multiwavelength 
observations and analysis of this bright afterglow and similar events (Harrison et al. 2001, 
Panaitescu & Kumar 2001, 2002, 2003, Willingale et al. 2004) established that afterglows 
were broadly consistent with the expected synchrotron spectrum and temporal evolution. 

If the relativistic outflow is collimated in the form of a jet then we expect to see an 
achromatic break in the decay at time tj days after the burst when the edge of the jet becomes 
visible (Rhoads 1997, 1999). Many optical observations of GRB afterglow decays exhibit a 
break a few days after the initial burst, which is identified with a jet break, consistent with 
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a collimation angle, 9j ~ 3 — 40 degrees (Frail et al. 2001, Bloom et al. 2003). Assuming 
the fireball emits a fraction rj^ of its kinetic energy in the prompt 7-ray emission and the 
circumburst medium has constant number density n the collimation angle is given by 



where z is the redshift, and E iso is the total energy in 7-rays in units of 10 52 ergs calculated 
assuming the emission is isotropic (Sari et al. 1999). The collimation-corrected energy is 
then E 1 = Ei SO (l — cosdj) and this shows a tight correlation with the peak energy of the 
spectrum in the source-frame, E^ r e c ak oc E®- 7 (the Ghirlanda relation: Ghirlanda et al. 2004). 
Jet breaks seen in the optical should also be observed, simultaneously, in the X-ray band. 

Prior to the launch of Swift (Gehrels et al. 2004, Burrows et al. 2005) both X-ray 
and optical follow-up observations of GRBs and their afterglows were limited to late times 
greater than several hours and often a day or more after the GRB trigger. Since launch, 
Swift has detected an average of 2 GRBs per week and we now have a sample of over 100 
GRBs for which we have quasi-continuous coverage in the X-ray band in the range ~ 100 
to ~ 10 6 seconds after the initial trigger. The aim of this paper is to compare the observed 
X-ray afterglows with the expectations of the standard model. One possible approach is 
to correlate the behaviour seen in the X-ray band with simultaneous optical measurements. 
Panaitescu et al. (2006) present an analysis for 6 GRBs detected by Swift noting that, 
contrary to expectation, temporal breaks in the X-ray band were not seen simultaneously in 
the optical. Another approach is to use the Ghirlanda relation to predict the time of expected 
jet breaks for Swift GRBs for which we have redshifts and then look to see if such breaks 
were observed. Sato et al. (2006) applied this to 3 bursts without success and concluded 
these bursts indicate a large scatter in the Ghirlanda relation. A third possible approach, 
presented here, is to make a systematic statistical study of the structure and evolution of a 
large sample of X-ray decay curves (including the method employed by Sato). 



An analysis of a sample of 40 X-ray decays observed by Swift by O'Brien et al. (2006), 
demonstrated that they all followed a similar pattern comprising an exponential decay in 
the prompt phase which relaxes to a power law decay at a time T p . In most cases this initial 
power law decay flattens into a plateau or shallow decay which then gradually steepens 
and establishes a final afterglow power law decay at time T a . Fig. [1] shows a schematic of 
the decay profile and the disposition of T p and T a . Such behaviour is consistent with the 
presence of two emission components that overlap in time; a short duration prompt emission 




(1) 



2. The functional form of X-ray decays seen by Swift 
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followed by an initial power law decay and designated by the subscript "p" and a much longer 
duration low luminousity afterglow component which starts as a slowly decaying plateau and 
ends with a steeper powerlaw, designated by the subscript "a". The analysis reported by 
O'Brien et al. (2006) concentrated on the properties of the prompt "p" component and 
produced an estimate of T p using a scaled version of each X-ray light curve. In this paper 
we turn our attention to the later development of X-ray light curves and employ a function 
fitting procedure to estimate the parameters associated with both the prompt and afterglow 
components. 

We have found that both components are well fitted by the same functional form: 
f c (t) = F c exp (a c - exp (^) , t < T c 



f c (t) = F c ( -J exp , t > T c (2) 

The transition from the exponential to the power law occurs at the point (T c , F c ) where the 
two functional sections have the same value and gradient. The parameter a c determines both 
the time constant of the exponential decay, T c j a c , and the temporal decay index of the power 
law. The time t c marks the the initial rise and the maximum flux occurs at t = \J t c T c / a c . 

Having established this generic behaviour we have fitted the X-ray decay curves of all 
107 GRBs detected by both the BAT and XRT on Swift up to August 1st 2006 using two 
components of the form f(t) = f p (t) + f a {t). Parameters with suffix p (T p ,. . .) refer to 
the prompt component and those with suffix a (T a ,. . . ) the afterglow component. Fig. [1] 
illustrates the functional form of the two components. 

The X-ray light curves were formed from the combination of BAT and XRT data as 
described by O'Brien et al. (2006). The conventional prompt emission, seen predominantly 
by the BAT, occurs for t < T p and the plateau/shallow decay phase, seen by the XRT, at 
t < T a . The fits were produced in two stages. The first stage used the BAT trigger time as 
time zero, to- in this fit the term exp(— t p /t) was included in the prompt function f p so that 
a peak position was found for the prompt emission. This peak time was then used as time 
zero and a second fit done with t p = (i.e. without an initial rise in the prompt component). 
Following this two stage procedure ensures that the prompt power law index fitted, a p , is 
referenced with respect to the estimated peak time rather than the somewhat arbitary BAT 
trigger time. In most cases the time of the initial rise, t a , of the afterglow component, f a {t), 
was fixed at the transition time of the prompt emission t a = T p . In a few cases this was 
shifted to later times because a small dip was apparent in the light curves before the start 
of the plateau or the plateau started particularly early. There was no case in which the 
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two components were sufficiently well separated such that this time could be fitted as a free 
parameter, i.e. we are unable to see the rise of the afterglow component because the prompt 
component always dominates/persists at early times and t a could be much less than T p for 
most GRBs. Many of the decays exhibit flares towards the end of the prompt phase, during 
the initial power law decay, on the plateau and even in the final decay phase. All large flares 
were masked out of the fitting procedure. Although apparently bright, such flares account 
for only ~ 10% of the total fluence in most cases. 

Chi-squared fitting was performed in log(flux) vs. log(time) space using the parameters 
log 10 (T p ), log 10 (T o ) and logs of the products, i.e. log 10 (F p T p ) and log 10 (T a T a ). The error 
estimation therefore produced a statistical error on the product of flux and time directly and 
these products could then be used to calculate the fluence and an associated fluence error 
in each of the components. The fluences of the prompt exponential and prompt powerlaw 
decay phases are 

F T 

fkxp = -^(exp(a p ) - 1) (3) 

Olp 

J- \ \ ^max J J 

where t max (> T p ) is the end of the light curve or some late time when the decay is deemed to 
have terminated. If a p >> 1 then t max can be set to infinity. The fluence of the exponential 
phase in the afterglow component is reduced by the initial exponential rise factor exp(— t a /t) 
and is given approximately by 

FT/ { t 



fhxp = [exp(a a y 1 - — ) ) - 1 ) (5) 



The inclusion of the exponential rise term has negligible effect on the fluence of the decay 
phase. Another way of viewing a p (or a a ) is that it controls the ratio of fluences seen from 
the exponential phase, t < T p (or T a ), and the decay phase, t > T p (or T a ). If the peak time, 
t p , is zero and t max — > oo then the ratio of the fluences for the prompt component is 

fhxp _ (exp(q p ) - l)(q p - 1) 

f^dec &p 

and flexp/ j I dec — 1 when a p = ot\ = 1.446. If a p « a\ then the decay is slow and most of 
the energy appears for t > T p in the power law decay. If a p » a>i then the decay is fast 
and most of the energy appears for t < T p in the early exponential phase. Fig. [1] shows the 
fluence ratio as a function of a p . A very similar expression holds for the fluence ratio of the 
afterglow component but this includes a minor adjustment because of the initial rise in the 
exponential phase, t <t a . 
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Table [T] lists the fitted parameters for all the GRBs in the sample. The type of decay 
fit (D) is also listed in Table [31 Of the 107 decays fitted, 85 required two components in 
which the afterglow component was dominant at the end of the observed light curve (D=l). 
In these cases the decay index of the prompt phase, a p was usually greater than the a a the 
final decay index of the afterglow. In a further 6 cases two components were required but 
the second component appeared as a hump in the middle of the initial decay and the prompt 
component reappeared and dominated again towards the end (D=2); these are marked* in 
Tabled) For these objects the prompt decay is slow, a p < 1.5, and it is usually the case that 
a p < a a . The remaining 16 required only one component (D=3) and did not exhibit a plateau 
phase with a subsequent power law decay. In 9 of these cases this was probably because the 
object was faint and the prompt emission faded below the XRT detection threshold before a 
putative plateau could be recognised. Fig. [2] shows typical examples of all these fitted types. 

For 99 of the 107 GRBs the latter stages of the light curve are well represented by the 
one or two component functional fit described above. For those with two components the 
plateau gradually steepens in the exponential phase, t <T a , and relaxes into a simple power 
law for t > T a . For those with one component the prompt emission turns over into a final 
power law decay at t > T p . However, in 8 cases there is clear evidence for a late temporal 
break. For these objects two extra parameters were included in the fit, a final break at 
time Tf, and a decay index a& for t > T&. Examples of these are also shown in Fig. [2] and 
Table [2] lists the fitted parameters. Such a late break usually occurs in the final decay of the 
afterglow component but in 3 GRBs, 060105, 060313 and 060607A, the late break is seen 
in the power law decay of the prompt component. For 060607A the break in the prompt 
component may have occured much earlier and the final break could be coincident with the 
end of the plateau, T a (see Fig. However, there is no doubt that a break occurs near the 
end of the light curve and the decay after this break is very steep. 

Fig. [3] shows the distribution of T a vs. T p for those objects with two component fits. 
There is no correlation between these timesJl] The frequency distributions of these times are 
shown in the bottom left panel. The same figure shows the distribution of afterglow fluence 
(the total fluence from the afterglow component) vs. the prompt fluence calculated as the 
sum of the T 90 fluence seen by the BAT and the fluence from the initial decay fld ec calculated 
from the XRT flux using the equation above. The dotted line indicates those objects for 
which the afterglow fluence is equal to the prompt fluence. There are a few objects on or 
just above this line while the rest are well below. There is a general trend that high prompt 
fluence leads to high afterglow fluence, as might be expected, but the scatter about this 



1 The error bars plotted in Fig. [3J and all subsequent plots involving fitted values, are 90% confidence 
ranges. 
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trend is large. This confirms the result from our earlier analysis (O'Brien et al. 2006) but 
for a larger sample. The frequency distributions of the fluences are shown in the lower right 
panel of Fig. [3J 

3. Spectral evolution 

Spectral fitting with XSPEC (Arnaud 1996) version 11.3.2 was used to determine the 
spectral index in the prompt phase (f3 p , from the BAT data), the prompt decay (/3 p d from 
the XRT data), on the plateau (p a , XRT data) and in the final decay (Pad, XRT data) for 
t > T a . In some cases the coverage was poor and/or the count rate low so it was not possible 
to separate (3 p d and (3 a or (3 a and /3 a d- For the weakest bursts it was only possible to derive 
(5 P and one spectral index from the XRT, (3 p d- When fitting late time spectra the absorption 
was fixed to the early time fitted values (both Galactic and intrinsic components) so that 
errors on the late time spectral indices were minimised. Table E] lists these spectral indices 
for all the GRBs in the sample. The ± ranges quoted are at 90% confidence. The decay fit 
type D is also listed. When D=3 there is no 2nd component and hence no plateau. However, 
for some of these afterglows late time XRT data are available and a late spectral index could 
be derived independently from (3 pd . These late spectral indicies are listed in the f3 a d column. 

Fig. H] shows the distribution of the afterglow plateau spectral index f3 a vs. the prompt 
(BAT over T90) spectral index j3 p . The range of indices from the prompt emission is large, 
—1.0 to 2.2 while the afterglow range is smaller, 0.4 to 2.3. Those objects for which the 
prompt emission is especially soft ((3 P > 1.8) or hard (f3 p < 0.1) evolve to produce an 
afterglow in the narrow range 0.7 < (3 a < 1.2. The frequency distributions of all the spectral 
indices are shown in bottom panels of Fig. |H The same figure shows the distribution of a a 
vs. a p for decays D — 1. Again, there is no correlation but the range of decay indices for the 
prompt component is large, 1.0 to 6.0, while the range for the 2nd afterglow component is 
much smaller, 0.5 to 2.0 with one object at a a = 2.7. Note, there are 4 decays with a a > 3 
in Tabled] but these are all type D = 2 marked*. In the majority of objects a a < a p . 

4. The expected coupling between a and (3 in the afterglow decay 

The standard fireball model predicts that there should be a simple coupling between the 
temporal decay index in the final afterglow, a a , and the spectral index in the final afterglow, 
I3 a d (Sari, Piran & Narayan 1998). The exact form of the relationship depends on the density 
profile of the surrounding circum-stellar medium (CSM) and whether the decay is observed 
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before or after a jet break (e.g. Panaitescu & Kumar 2002); because the spectrum is being 
red-shifted as the jet is slowed down by the CSM and after a jet break the peak flux of 
the observed synchrotron spectrum is also decaying with time. The top left panel of Fig. [5] 
shows the temporal decay index of the afterglow component, a a , vs. f3 a d the spectral index 
in the final decay for those objects where t max > T a and for which we have a significant 
measurement of (3 a d at t > T a . There is no correlation between a and j3 and if there is any 
trend at all it is that some of the faster decays (large a) occur for the smaller (3 values. The 
lower band shown indicates the region of a pre-jet break coupling predicted by the model and 
the upper band the region expected post-jet break. In each case the lower edge of the band 
corresponds to the model in which the X-ray frequency is below the cooling break frequency 
and the CSM density is constant. The upper edge corresponds to the X-ray frequency above 
the cooling break frequency and a wind density oc r~ 2 . Of the 70 objects plotted, 36 lie 
below the expectations of the standard model in the bottom right of the plot. For 17 the 
upper limit of the 90% confidence region in (3 — a doesn't intersect the pre-jet break band. So 
for 50% of GRBs the spectral index of the afterglow is too large to produce the observed 
temporal decay index. This can only occur under the model if there is significant energy 
injection such that the peak of the spectrum (or normalisation) is boosted fast enough to 
counteract the drop expected from the change in red-shift as the outflow slows down. It 
was pointed out by Nousek et al. (2006) that the spectrum and decay of the plateau phase 
of several Swift GRBs was inconsistent with the expectations of the model and that energy 
injection during this phase was a possible explanation. The current analysis shows that the 
same is also true for many objects during the subsequent decay phase, after the plateau. The 
lower left-hand panel of Fig. [5]shows a a — (3/2)/3 ad plotted as a function of T a . This function 
will be zero for pre-jet break afterglows under the standard model (with uniform CSM and 
the X-ray frequency below the cooling frequency) and negative for afterglows with a value of 
a a which is too small compared with the /3 ot j. The horizontal dashed line is the pre-jet break 
expectation if the X-ray frequency is above the cooling frequency. Consideration of GRB 
efficiencies (Zhang et al. 2006) indicates that more than 60% of afterglows in the sample 
described by O'Brien et al (2006) do lie above the cooling break. If this is the case for the 
present, larger, sample a significant fraction of the GRBs plotted as dots (pre-jet break) are 
also outside expectation. The figure also indicates that there is no correlation of this function 
with the time at which the final decay starts, so if energy injection is the explanation it is 
occurring at large times, > 10 5 seconds as well as much earlier times. Of the remaining 
~ 50% of objects plotted on the top left-panel of Fig. [5j 26 lie within the pre-jet break 
band while 8 lie above this close to or within the post-jet break region. It is conceivable 
that some of the anomalously large (3 values could be associated with the fitted absorption 
(both Galactic and intrinsic). We checked this possibility but found no correlation between 
the position of an afterglow on the (3 — a plane and the fitted Nh- 
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The top right-hand panel of Fig. [5] shows the change in spectral index into the final 
decay from the plateau, A/5 = f3 ad — (3 a plotted vs. the plateau spectral index (3 a . Here there 
is a trend. If (3 a is small then A/5 > and the afterglow gets softer in the final decay. If f3 a 
is large then A/5 < and the final decay is harder. The scatter of index f3 a d is smaller in 
the final decay, being confined to the range 0.0 to 2.2 as indicated on the top left plot. The 
gradual narrowing of the spectral index range as the afterglows develop can also be seen in 
the frequency distributions shown in the lower panels of Fig. HJ One object, GRB060218, 
has an anomously large final spectral index, j3 ad = 3, but this GRB was very peculiar in 
many respects, in particular for having a significant thermal component in the early X-ray 
spectrum, Campana et al. (2006). No (3 a was derived for this afterglow because the plateau 
is largely obscured by unusual, persistent, prompt emission. The trend in the change in 
spectral index from the plateau into the power law decay is independent of the position in 
the P — a plane occupied by the final decay. 

The above discussion has considered those objects with a 2nd afterglow component 
that dominates in the later regions of the X-ray light curve (D=l). 22 of the X-ray decays 
required no 2nd component (D=3) or have a 2nd component which produced a hump in 
the decay but faded towards the end (D=2). In these cases the initial decay of the prompt 
component dominates at late times. The bottom right-hand panel of Fig. [5] shows the 
correlation between the late temporal decay index (a p ) and spectral index of these objects. 
The behaviour is very similar to objects (D=l) shown in the top left panel. There is no 
obvious correlation and 8 of the 21 object fall below the expected correlation in the bottom 
right. Also plotted in this panel are the f3 a d,(Xb values for the late breaks listed in Table 
H Five of these lie within the pre-jet break band but three, GRB050814, GRB060105 and 
GRB060607A, are within the post-jet break region although the errors on the decay index 
after the late break temporal break, o^, are rather large. The complete X-ray light curve for 
GRB050814 containing the late temporal break is shown in the upper right panel of Fig. [5J 
Further analysis of GRB060105 is provided by Godet et al. (2006a). GRB06067A is similar 
to GRB060105. 



5. Isotropic energy of the prompt and afterglow components 

Fig. E] shows the correlation between the fluences in the prompt component (including 
the BAT Tg 15 to 150 keV and the initial decay in the XRT 0.3 to 10 keV) and the afterglow 
component. For those bursts for which we have a redshift (z values listed in Table |3]) these 
fluences can be used to estimate the equivalent isotropic energy. We assumed a cosmology 
with H Q = 71 km s _1 Mpc" 1 , A = 0.27 and Q = 0.73 and calculated E iso over an energy 
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band of 1 to 10000 keV in the rest frame applying a k-band correction independently for 
the BAT and XRT data. For many objects we don't have a measured value of the peak 
energy in the spectrum. In such cases we assumed E peak = 116 keV (the median value for 
Swift bursts) and a spectral index of (3 = 2.3 for E > E pea ^. The mean E pea k for BATSE 
(pre-Swift) bursts was 235 keV (170-340 keV 90% range) (Kaneko et al. 2006) which is 
significantly larger than for Swift bursts, mean 138 keV, median 116 keV. However the mean 
redshift for Swift bursts considered here is 2.46 while for pre-Swift bursts it is ~ 1.52 so we 
expect the mean/median E pea k for Swift bursts to be ~ 0.6 the BATSE value. The mean 
value for the spectral index above E pea k for BATSE bursts was 2.3 (Kaneko et al. 2006). 
Fig. [6] shows the equivalent isotropic energy for the two components. The symbols used are 
the same as in Fig. [5j It is clear that the total energy in the afterglow is not correlated 
with the position of the final afterglow in the (3-a plane. Note that the correlation between 
the equivalent isotropic energy components evident in Fig. [6] is real but not particularly 
significant since it arises from applying the measured redshift in both axes. In many objects 
the total energy seen in X-rays from the afterglow is significant compared with the 7 and 
X-ray energy seen from the prompt component. The dotted histogram in the right panel 
shows the distribution of E iso values for GRBs observed by instruments pre-Swift taken from 
the tabulations in Frail et al. (2001), Bloom et al. (2003) and Ghirlanda et al. (2004). The 
maximum isotropic energy in the Swift sample is similar to the maximum seen previously, 
« 10 54 ergs, but the distribution of energies seen by Swift is broader, has a lower mean and 
extends to a lower limit of ~ 10 47 ergs. 



6. Temporal breaks 

The visibility of the 2nd component used here to fit the X-ray decay curves depends on 
the relative brightness of the prompt emission decay compared with the afterglow plateau 
and the times T p and T a . If T a is not long after T p then the end of the plateau phase is not 
visible, as is the case for the decay shown at the top left of Fig. [2j However, for 64 of the 
91 GRBs which required 2 components in the fit the end of the plateau is visible, as is the 
case for the exemplar GRB shown at the top right of Fig. [2j In such cases the plateau gets 
slowly steeper towards the end of the exponential phase and eventually relaxes to a power 
law. There is often no definitive or sharp break but the time T a is a robust measure of where 
this transition occurs, taking into account all the data available. Thus, the fitting provides 
91 afterglow break times, T a , from a total of 107 objects. 

In some cases it may be that any jet break time associated with the edge of a putative jet 
becoming visible occurs at or before T a . In such decays we expect the subsequent afterglow 



- 11 - 



to lie somewhere in the top left of the (3 — a plane shown in Fig. [5] and the 8 candidates for 
such cases are shown as star symbols on this figure. Note that the error bars shown in Fig. 
[5] are at 90% confidence and there is one object lying in the lower pre-jet break band which 
is also consistent with the upper post-jet break band. The time T a for these GRBs is not a 
pure jet break time since in all cases the T a marks the end of the plateau phase which does 
not behave as a pre-jet afterglow (e.g. Nousek et al. 2006). The only evidence for a jet break 
having occurred in these 8 candidates is that the a and (3 values of the subsequent decay 
have the right relationship for a post-jet break afterglow. Since 50% of afterglows don't 
agree with the expected alpha-beta value, including the pre-jet break values, the argument 
in favour of jet breaks in these cases is weak. However, T a is a reasonable estimate of where 
any jet break may have occurred. Details of these 8 potential jet breaks are given at the top 
of Table H 

As discussed above, an additional 8 GRBs required a late break to fit the data as 
illustrated at the bottom of Fig. [21 The positions of the final afterglows in the (3 — a plane 
for these GRBs are shown in the bottom right of Fig. [5] (plotted as triangles). Only 4 of 
these objects are good candidates for jet breaks and these are listed at the bottom of Table 
H] and illustrated in Fig. [HJ The jet break time ranges plotted were estimated using the 
Ghirlanda relation assuming E peak = 116 ± 50 keV (there are no measured values of E peak 
for these bursts). For GRB060607A, in the bottom right panel, the late break occurs slighly 
earlier the allowed band indicating that, for this break to be consistent with Ghirlanda, the 
E peak should be < 66 keV. So, in summary, out of 72 afterglow breaks identified by the 
fitting procedure only 12 are followed by an afterglow which is consistent with post-jet break 
conditions and only 4 of these are isolated breaks independent of the end of the plateau 
phase. The remaining 60 have slow X-ray decays (low a) and/or soft X-ray spectra (high 
0). 

The interest in jet breaks stems from Equation [T] which can provide an estimate of 
the jet angle 9j and hence the collimation-corrected energy E T For calculation of all the 
subsequent E 1 values considered below we have assumed = 0.6 (or equivalently n = 3 
cm -3 and r? 7 = 0.2). r) 7 = 0.2 has been widely assumed in pre-Swift analysis although recent 
work using Swift data (Zhang et al. 2006) indicates that the efficiency can be determined 
with more accuracy. However, since E 1 oc (m^) 1 / 4 the resulting collimated energy is fairly 
insensitive to these parameters. For 9 of the objects in Table H] we have redshifts and can 
estimate E 1 using the jet break time. These are shown in Fig. [7J The symbols are the same 
as in Fig. [5J Also shown are the values derived for pie-Swift GRBs using the tabulations in 
Frail et al. (2001), Bloom et al. (2003) and Ghirlanda et al. (2004). For 31 others we have 
redshifts (as listed in Table [3]) but there is no break between T a and the last data point at 
tmax- For these decays we can calculate a range of E 1 which is excluded by the decay. These 
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ranges are plotted in the top right-hand panel of Fig. [7] with the energy corresponding to 
T a at the lower end of each range and to t max at the upper end. For many afterglows the 
excluded range of E y covers a substantial fraction of the pre- Swift distribution. The lower 
panels of the same figure shows the respective frequency distributions. The 3 objects with 
measured late breaks are in good agreement with the peak of the pre- Swift E 1 distribution 
published by Frail and Bloom. Of the remaining 6 objects in which T a has been identified 
with a possible jet break 2 lie within the lower wing of the pie-Swift distribution and 4 
lie below. The distribution of E n values derived from the remaining T a values is similar in 
shape to the pie-Swift distribution calculated using optically observed jet break times but 
is offset to lower energies by a factor of ~ 34. This corresponds to an average jet break 
time which is a factor « 110 smaller or a jet angle which is a factor w 10.5 smaller. The 
times T a derived from the X-ray decay curves are, on average, a factor of « 80 smaller than 
the optical jet break times observed for pie-Swift GRBs. The peak of the distribution of 
E 1 values calculated from t max is significantly higher than the pie-Swift distribution because 
many of the X-ray decays extend to later times without a temporal break. 

T a is the only measure we have (or require) to specify the timescale of the 2nd afterglow 
component. It is clear from the discussion above that, in almost all cases, this does not 
represent a jet break time. However, given that it is the only time extracted from the X-ray 
data it may be related to the previously observed optical jet break times. There are 14 
GRBs for which we have a T a value, a redshift and a measurement of the peak energy of 
the 7-ray spectrum, E pea ^. For these objects we can construct a Ghirlanda type correlation 
between the peak energy in the rest frame E pea k{z + 1) and the collimated beam energy 
calculated assuming tj = T a . This is shown alongside the conventional Ghirlanda relation in 
Fig. [91 The X-ray measurements show a similar behaviour, not as statistically significant as 
the optical version but with about the same gradient and an offset in E 1 by a factor of ~ 26. 
For these cases the X-ray break times, T a , are a factor of ~ 90 less than pie-Swift optical 
break times. Unfortunately none of the GRBs with potential X-ray jet breaks identified and 
plotted in Fig. [8] has a measured E pea k value so these can't be included on Fig. [9j However, 
we can calculate lower limits to E 1 , assuming tj = t max , for those decays without a late 
break but with an E pea k and redshift measurement. These are included on Fig. [9J They all 
lie to the right of the Ghirlanda correlation. There is currently one afterglow observed by 
Swift, GRB050820A, for which we have derived a time T a and the optical coverage extends to 
~ 100 x T a (Cenko et al. 2006). In this case the optical decay does appear to have a break at 
about the expected time, tj ~ 100 x T a . The Ghirlanda correlation was originally formulated 
by considering the jet model and the collimated jet energy however, the correlation derived 
here, using T a from the X-ray afterglows, has more in common with the model independent 
analysis described by Liang & Zhang (2005). The Liang-Zhang E iso — E pea k—tb y0 pt correlation 
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may well be related to the E iso — E pea k — T a correlation plotted in Fig. [9j 

As a natural consequence of the above discussion, and following Sato et al. (2006), we 
can adopt a different approach and look for X-ray decays in which a predicted jet break 
is absent. Apart from GRB050401 (or XRF050401, Mangano et al. 2006, Sakamoto et al. 
2006), GRB050416A and GRB050525A already considered by Sato et al. (2006) there are 5 
other X-ray decays monitored by Swift which come under this category. Light curves for all 
8 are shown in Fig. [TDJ In each case the expected range for the jet break time, tj, is shown 
as the shaded area. Each decay was followed for at least a factor of 100 times longer than 
T a without any indication of a break in the temporal decay index. For GRB060206 there is 
some indication that the X-ray light curve is flattening for t > 10 6 seconds. This is probably 
due to systematic errors in the background subtraction when the afterglow is very faint 
and/or contamination by a faint, nearby, background object. GRB050525A, GRB050820A, 
GRB060510B and GRB060729 are cases where there is a large flare at the end of the prompt 
phase and the late BAT and/or early XRT data are not well fitted by the prompt function. 
However, the plateau and subsequent power law decay are well represented by the afterglow 
component in all cases. The late afterglow of GRB060729 has currently been monitored for 
77 days without any indication of a break. 

Optical data are available for 5 of the final power law decays of these afterglows; 
GRB050525A (Blustin et al. 2006), GRB050820A (Asian et al. 2006, Cenko et al. 2006), 
GRB060206 (Stanek et al. 2006, Monfardini et al. 2006), GRB060707 (de Ugarte Postigo 
et al. 2006, Jakobsson et al. 2006) and GRB060729 (Grupe et al. 2006). Optical jet breaks 
are not seen in any of these afterglows although the decays of GRB050525A, GRB060206 
do gradually steepen at late times and there is a one final late measurement at ~ 3 x 10 6 
seconds for GRB050820A which indicates that the optical decay has turned down (HST data 
presented in Cenko et al. 2006 as already mentioned above). We have fitted these optical 
afterglows with a simple power law over the period contemporary with the final X-ray decay 
(t > T a ) and the results are shown in Table [5j For 4 afterglows the X-ray and optical decay 
indices are consistent. For GRB050820A the optical index is formally significantly lower than 
the X-ray index but there is structure in the optical decay curve which is not well modelled 
by a single power law fit. For all 5 the X-ray and optical decays are remarkably similar and 
shallow with decay indices in the range 0.75-1.5 and they all have a, (3 reasonably consistent 
with a pre-jet break afterglow. 

Two other objects, GRB050822 (Godet et al. 2006b) and GRB060218 (Campana et al. 
2006) also have Swift X-ray light curves with extended coverage of a long power law decay 
without a break. Both these objects have soft spectra and are more reasonably considered 
as XRFs with E peak < 25 keV. For GRB060218 we have a measured redshift of 0.0333 
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(Mirabal & Halpern 2006). Adopting a range of E peak values consistent with an XRF (and a 
range of redshift values for GRB050822) it can be shown that a jet break should have been 
seen at some point in the X-ray light-curves of both these objects but none was observed 
(Godet et al. 2006b). For GRB060218 radio observations have shown that the fireball was 
probably isotropic (Soderberg et al. 2006) and similarly, in the cases of GRB050416A and 
GRB050822, it is possible that the jet opening angle was wide (> 20 degrees). For the XRFs 
a large opening angle and consequent absence of a jet break is not unreasonable because 
E iso is typically rather low. There are also 9 further GRBs (050318, 050603, 050802, 060108, 
060115, 060418, 060502A, 060614 and 060714) for which we can predict a range of jet break 
times and no break is seen. For these objects the lowest predicted tj is well before the last 
observed time but the upper limit on tj is after the last observed time so we cannot rule out 
the possibility that a break occurred too late (and too faint) to be seen by the Swift XRT. 
There is a suggestion of a late break for GRB060614 but inclusion of such a break is not 
statistically significant under the current analysis. 

An important property of a jet break is that it should be achromatic, occurring across 
the spectrum with the same temporal profile. Panaitescu et al. (2006) show that temporal 
breaks for 6 GRBs seen by the Swift XRT are not present in optical data which span the 
same period of time. For 5 of these, GRB050802, GRB050922C, GRB050319, GRB050607 
and GRB050713A the break in the X-ray decay is fitted as T a in the analysis presented 
above. Therefore for these objects optical data do not follow the X-ray profile modelled by 
the 2nd afterglow component and the observed X-ray break is unlikely to be a jet break. The 
remaining object, GRB050401, is fitted using 2 components but the 2nd component appears 
as a hump in the X-ray decay and the 1st component has a slow decay, a p = 1.12 ± 0.05, 
that dominates near the end. So in this case the optical data seem to follow the behaviour 
of the power law decay of the 1st prompt component and not the 2nd component. Similar 
behaviour is exhibited by GRB060210 which has the same relationship between the 1st and 
2nd X-ray component, a p = 1.0 ± 0.07, and the available optical data (Stanek et al. 2006) 
follow the decay of the 1st component. Conversely, there are examples of GRB decays for 
which the X-ray and optical profiles do follow a similar pattern. Optical data for GRB060206 
(Stanek et al. 2006) have a profile which closely follows the 2nd component of the X-rays. 
Stanek et al. argue that for this GRB there is a break which can be seen simultaneously in 
optical and X-rays but the present analysis finds no such X-ray break. Instead there is a 
gentle curvature in the X-ray (and optical) decay which is modelled by the slow transition 
from exponential to power law in the profile of the 2nd afterglow component. GRB050525A 
is another example in which curvature in the later stages of the X-ray and optical can be 
modelled as a break (Blustin et al. 2006). The present analysis finds no such break in the 
later stages of the X-ray decay of this burst either. 
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7. Conclusions 

The X-ray decay curves of 107 GRBs observed by Swift have been fitted in a systematic 
way using the simple functional form given by Equation [2j They all require a prompt 
component with parameters T p , F p and a p and 85 require a 2nd afterglow component with 
parameters T a , F a and a a . The parameters and associated confidence limits are all listed 
in Table [H T p is similar to the familiar Tgo burst duration as discussed by O'Brien et al. 
(2006). a p is a prompt parameter which was unavailable before the Swift era and indicates 
how fast the prompt emission is decaying as also discussed by O'Brien et al. (2006). The 
product T p F p combined with a p is a measure of the prompt fluence (see Equations [3J and @J . 
a p also determines the distribution of energy. If a p < 1.446 then more energy is emitted for 
t > T p during the prompt power law decay phase and if a p > 1.446 more energy is emitted 
during the prompt exponential phase t < T p . T a is the time when the final afterglow power 
law decay starts, a a is the index of this final decay and also controls the curvature of the 
proceeding plateau phase and the product T a F a along with a a combine to give a fluence 
of the afterglow component (sum of Equations [3J and H] replacing subscripts p with a). a a 
controls the distribution of energy in the 2nd functional component in the same way as a p 
controls this distribution in the prompt component, as described above. In addition, spectral 
fitting has provided 7 and X-ray spectral indices over the phases of these functional fits as 
listed in Table [31 The combination of Tables [Hand [3J provides a rich data base for comparison 
with theoretical models. 

The functional form used for the fits (Equation [2]) is empirical rather than based on 
a physical model and it doesn't accomodate flares. However, the exponential rise at time 
t c and transition from exponential to power law at T c are reminiscent of the theoretical 
discussion of the development of an afterglow given by Sari (1997). Employing exponentials, 
rather than power law sections with breaks, provides a curvature which neatly fits the data 
and produces a remarkably good representation of the underlying X-ray decay profile with 
the minimum number of parameters. In any physical model the number of parameters that 
could influence the shape of the X-ray light curves is large but some combination or subset 
of these are likely to be represented by the values tabulated here. 

Most (?^ 80%) GRBs have a 2nd afterglow component which dominates at later times 
and this is probably the expected emission from the external shock. For these objects the 
prompt and afterglow components appear to be physically distinct. Bursts marked* in Table 
[1] (or D=2 in Table [3j have 2 component fits but the 1st prompt component dominates at 
later times. For these GRBs and those which only require 1 component in the fit (D=3 in 
Table [3j it is not clear where such an early prompt decay component comes from. It could 
be the external shock but if this is the case then the external shock is developing very early 
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(t < T p ) and the prompt decay emission from this shock is distinct from the more common 
external shock emission seen in the 2nd component. In GRBs marked* (D=2 in Table [3]) 
both these external shock components are seen. For these objects it is not so obvious that 
the prompt and afterglow components fitted map directly to two physical entities. The 
hump in the decay curve could be a signature of the presence of different outflow structures 
(e.g. Eichler & Granot 2006) or a universal structured jet (e.g. Meszaros, Rees & Wijers 
1998). Two of the GRBs, GRB051210 and GRB051221B, fitted by just one component (D=3 
in Table [3]) could be "naked GRBs" for which there is no external shock and the prompt 
decay arises solely from high latitude emission (Kumar & Panaitescu 2000). Such GRBs are 
expected to have a p - (3 pd = 2. This difference is 1.81 ± 0.2 for GRB051210 and 1.73 ± 1.0 
for GRB051221B consistent with a high latitude interpretation. Neither of these bursts were 
detected after the 1st orbit post trigger. GRB051210 has an upper limit from the 2nd orbit 
consistent with the steep power law decay. For GRB051221B upper limits were obtained 
from the 2nd orbit and much later, t = 4.6 x 10 5 seconds, both again consistent with a steep 
prompt decay and no afterglow or plateau. The remaining 14 single component prompt 
decays were not steep with mean a p = 1.35 ± 0.22 and significantly lower a p — (3 p d values 
with a mean of 0.6 ± 0.4. Although a rise time, t a , was included in the afterglow function fit 
in most cases it was set at some arbitarily early time when any increase in this component 
would be obscured by the prompt emission. Thus, the functional fits provide very little 
evidence for any observational signature of the "onset" of the afterglow or a rising afterglow 
component predicted in some off-axis models (Eichler & Granot 2006). The afterglow could 
be established very early. Since the outflow is presumed to be moving close to the speed of 
light there may be very little time delay between prompt emission and emission from the 
external shock when the outflow begins to decelerate. 

There is no tight correlation between the prompt component parameters and the after- 
glow component parameters. There is a trend that the larger the prompt fluence the larger 
the afterglow fluence, as might be expected, and the afterglow energy is usually considerably 
less than the prompt energy. In some cases the afterglow energy is about equal to the prompt 
but we never see afterglows which are much more energetic than the prompt emission; see 
Figures [3] and El 

Almost all of the X-ray afterglows end up in the same region of the (3 — a plane, 
0.5 < (3 < 2 and 0.4 < a < 2, independent of the other parameters. However, « 50% of 
these afterglows are in the bottom right of the (3 — a plane with high (3 and low a, values 
not predicted by the standard fireball model. This may be due to persistent long lasting 
energy injection but this seems unlikely in such a high percentage of objects and at such late 
times. Furthermore the trend expected from the model, that high (3 values should give high 
a values and vice versa, is not seen at any level of significance. If this is because of energy 
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injection this energy injection is occurring in just the right number of afterglows and at just 
the right level such that the expected j3 — a coupling is hidden or masked. If, for example, 
energy injection were always present at some fixed level all afterglows would be shifted to 
lower a values but any correlation would be preserved. In fact, if energy injection occurred 
preferentially at late times post-jet break afterglows would be shifted to smaller a values 
but pre-jet break afterglows would not be altered and the overall correlation including both 
pre and post jet break afterglows might be tightened. The curvature or break which marks 
the end of the afterglow plateau phase at time T a is often accompanied by a small spectral 
change such that soft afterglows become somewhat harder and hard afterglows become softer 
bringing the final afterglow spectral index into the rather narrow range indicated above. 

For the 91 GRBs which require a 2nd afterglow component we derive time T a which 
marks the end of the plateau and the start of the final decay. For 64 cases we can see this 
transition/break in the data. Do we see any jet breaks? For 8 GRBs the final temporal decay 
and spectral index t > T a are consistent with the jet model after a jet break has occurred 
but for 6 of these objects for which we have a redshift the Ej values are low compared with 
the distribution derived from consideration of optical jet breaks and the implied jet angles 
are small. 8 X-ray decays have a late break, T& > T a , and 4 of these have final a and f3 values 
consistent with a post jet break afterglow. 3 of these have redshifts and E 1 values which 
are in accord with the pre-Swift distribution. On the other hand there are 11 decays with 
long-lasting coverage, t > T a , in which a jet break predicted using the Ghirlanda relation is 
definitely not seen and there are a further 9 objects in which a break is not seen but it may 
have occurred at a rather late time beyond the coverage provided by the sensitivity limit 
of the Swift XRT. 13 decays (in the top right and bottom left panels of Fig. [5]) have been 
identified as lying in the post jet-break region of the (3 — a plane (star symbols). For such 
decays the electron energy distribution index is expected to be equal to the decay index. 
Two of these decays, GRB060421 and GRB060526, have decay indices < 1.6 (see Table 0J, 
probably too low to be the same as the electron index. 

There are currently ~ 10 afterglows for which there are both X-ray and optical data 
available during the X-ray plateau and into the final decay (and there are a few more objects 
for which such data will soon be publically available). The X-ray and optical often follow 
a similar trend but temporal breaks which are seen in one band are not always seen in the 
other. This is probably because the end of the plateau and start of the final power law 
is marked by a continuous curvature such that the light curves in both bands are slowly 
getting steeper. Rather than fitting isolated temporal breaks independently in either band 
we suggest that simultaneous fitting of a profile as described by Equation [2] in an attempt to 
find a value for T a which is consistent with both the X-ray and optical data would be more 
illuminating. Such a joint fit would also provide a best estimate of the fluence ratio between 
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the X-ray and optical bands (or equivalently the X-ray/optical flux ratio). 

If we assume that T a is a jet break time tj and calculate E 1 for those objects with 
redshifts the distribution of E 1 which results is very similar to the distribution derived from 
optical jet breaks but is offset to lower energies by a factor pa 34 implying that any optical jet 
break should be seen at time tj ~ 110 x T a . Furthermore, for those GRBs for which we have 
an E pea k value for the 7-ray spectrum the peak energy in the rest frame is correlated with the 
collimated energy estimate E 7 , in the form of a Ghirlanda relation, but with the E 1 values 
offset to lower energy by a factor pa 26 which implies that tj pa 90 x T a . It appears that T a 
extracted from X-ray decays in the present analysis has properties related to tj derived from 
optical data. Whether this apparent connection is purely statistical in nature or has some 
deeper significance remains to be seen but it is doubtful that either T a or tj are actually "jet 
break" times. However, it is likely that the end of the plateau phase, T a , does depend on the 
total energy in the outflow, the collimation angle of the outflow and the density of the CSM 
and that the E iso — E pea k — T a correlation reported here is related to the E iso — E pea k — tb, op t 
relation discussed by Liang & Zhang (2005). 
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Fig. 1. — Left: The functional form of the decay and the fitted parameters. The prompt 
component (green) has no rise because time zero is set at the peak. The afterglow component 
(blue) rises at time t a as shown. Right: the ratio of the fluence in the exponential portion of 
the light curve to that in the power law decay as a function of a p . See Equations [3j H] and [6j 
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Fig. 2. — Examples of the fits to X-ray decay curves. The prompt and afterglow component 
functions are plotted in the same way as in Fig. [TJ Top panels: the most common type 
in which the 2nd, afterglow, component dominates at late times. Middle panel left: a 
two component fit in which the afterglow component forms a bump in the decay but the 
extrapolation of the prompt component decay dominates at late times. Middle panel right: 
a single component fit, requiring no afterglow component. Bottom panels: two examples of 
fits which include a late temporal break. Flares are plotted in red and were excluded from 
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Fig. 3. — Top left: The afterglow duration T a vs. the prompt duration T p plotted for all 
objects in which the 2nd component dominates the prompt component at late times. Top 
right: The afterglow fluence calculated by integrating the 2nd component in the XRT band 
0.3-10 keV vs. the prompt fluence calculated from the BAT T90 flux 15-150 keV plus the XRT 
flux 0.3-10 keV in the initial decay (fldec)- The dotted line indicates where the prompt and 
afterglow fluences are equal. Bottom left: Frequency distributions of T p (open histogram) 
and T a (shaded histogram). Bottom right: Frequency distributions of prompt fluence (open 
histogram) and afterglow fluence (shaded histogram). 
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Fig. 4. — Top left: Comparison of the spectral indices f3 a (2nd, afterglow component) with 
f3 p (prompt BAT). Top right: Comparison of the decay index for the 2nd component, a a with 
the prompt component, a p for decays type D — 1. The dotted lines indicate equality in both 
panels. Bottom left: Frequency distributions of the prompt component spectral indices, (3 P 
(dotted histogram) and (3 p d (shaded histogram). Bottom right: Frequency distributions of 
the afterglow component spectral indices, [3 a (dotted histogram) and f3 a d (shaded histogram). 
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Fig. 5. — Top left-hand panel: the temporal decay index of the afterglow plotted vs. the 
spectral decay index in the final afterglow decay. The lower band (yellow) indicates the 
area expected to be occupied by a pre-jet break afterglow and the upper band (orange) the 
area post-jet break. GRBs which fall in the pre-jet break region are plotted as circular dots 
(black), those which fall above this in the post-jet break region are plotted as stars (red) and 
those below the pre-jet break band are plotted as squares (green). The same symbols are 
used in the other three panels. Top right-hand panel: change in spectral index (3 a d — /3 a vs. 
the plateau spectral index (3 a . Bottom left-hand panel: a a — (3/2)/? ad vs. the duration of the 
afterglow T a . This function is expected to be zero for a pre-jet break afterglow in a uniform 
CSM with the X-ray frequency below the cooling break frequency. The dotted line indicates 
the zero if the X-ray frequency is above the cooling break frequency. Bottom right-hand 
panel: the temporal decay index vs. the spectral index for afterglow power law decay in 
objects which require only one component in the fitting or in which the prompt component 
dominates at the end of the X-ray light-curve. The j3 a d vs. for the late temporal breaks 
listed in Table [2] are also shown as triangles (blue) . The bands shown are the same as for 
the top left-hand panel. 
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Fig. 6. — Left: The isotropic afterglow energy vs. the isotropic prompt energy. Symbols 
indicate the position of the afterglow in the (3- a plane as in Fig. The dotted line indicates 
equality between the prompt and afterglow energies. Right: The distribution of isotropic 
energies. Filled histogram E prompt and solid line histogram E a f ter giow The dotted histogram 
is the distribution of Ei SO for GRBs observed before the Swift era as tabulated by Frail et al. 
(2001), Bloom et al. (2003) and Ghirlanda et al. (2004). 
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Fig. 7. — The jet energy E 1 calculated using equation [TJ Top left: The stars and triangles 
correspond to the jet break times listed in Table H] (symbols as Fig. The open circles 
are the values obtained using the tj and Ei SO tabulated by Frail et al. (2001), Bloom et al. 
(2003) and Ghirlanda et al. (2004). Top right: Ranges of E n excluded by the X-ray decays 
without breaks, t > T a . The lower points are from T a and upper from t max (the end of the 
light curve). Bottom left: The filled histogram shows the distribution of E 1 for the breaks 
listed in Table HI The dotted histogram is the Frail-Bloom- Ghirlanda sample. 
Bottom right: The dotted histogram is the Frail-Bloom-Ghirlanda sample. The solid line 
histogram is the distibution derived from T a and the dashed line histogram from t max . 
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Fig. 8. — X-ray decays in which there is potential jet break. The shaded area indicates 
the expected range of tj calculated using the Ghirlanda relation for GRBs with a measured 
redshift. None of these has a measured E pea k value so 116 ± 50 keV has been assumed (see 
Section EJ). The estimated E iso (10 52 ergs) and predicted tj (days) are listed on the left. 
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Fig. 9. — The open circles and solid line are the Ghirlanda relation. The points with 
error bars are Swift GRBs for which values of E 1 were derived using the Ghirlanda relation 
assuming tj = T a . The symbols plotted are the same as in Fig. [5] The dashed line is a fit to 
the Swift sample assuming the slope is fixed at 0.706, the best fit value from Ghirlanda et 
al. (2004). The lower limits of energy shown were derived using the latest observed time as 
tj for GRBs with lasting afterglow decays which show no break (the light curves plotted in 
Fig. [TO]) . -Ep ea fc values used here and for Figures [8] and [10] were taken from Cummings et al. 
(2005), Crew et al. (2005), Golenetskii et al. (2005, 2006) or they were derived from spectal 
fitting of the BAT data. 
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Fig. 10. — X-ray decays in which a predicted jet break is absent. E peak listed in keV, 
Ei SO in 10 52 ergs and tj in days. The shaded area indicates the expected range for tj. For 
GRB060510B and GRB060729 there is no measured E pea k value so 116 ± 50 keV has been 
assumed (see Section E]). 
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Table 1. Best fit parameters for the Swift X-ray decay curves together with upper and 
lower 90% confidence limits. Times are in seconds and fluxes in ergs cm -2 s _1 0.3-10 keV. 
t p is the prompt peak time with respect to the BAT trigger time used as time zero. *2 
component fits in which the 1st prompt component dominates towards the end of the 
X-ray light curve. +These GRBs have less than 50% of T 90 included in the fitting. 
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Table 2. Late temporal breaks with 90% confidence limits. Tf, is the break time in 
seconds. The decay index before the break was either a p or a a as shown depending on 
which component dominates at the end of the X-ray light curve. The final decay index 
after the late break is given by a&. Aa is the change in decay index across the break. For 
060607A the late break is close to the decaying section of the 2nd component so the Aa 

value maybe misleading/incorrect. 
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Table 3. Spectral indices; (3 P prompt phase, (5 p d initial decay, (3 a plateau and (3 a d final 
decay after the plateau. Limits are 90% confidence. Column D is the type of decay fit from 
Table [TJ 1 - components p and a required with component a dominant at end, 2 - 
components p and a required with component p dominant at end, 3 - only component p 

required. Column z is the measured redshift taken from 
http:/ /swift. gsfc.nasa.gov/docs/swift/archive/grb_table/. References for all the redshifts 

are provided on this WWW data table. 
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060421 
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0.49 ±0.05 


U.54 + U.33 


U.41 + U.39 


U.33 + U.37 




060427 


1 


0.90 ±0.17 


1.90 ±0.21 


U.54±U.34 


1.15 ±0.17 




060428A 


1 


1.04 ±0.07 


2.35 ±0.25 


1.20 ±0.15 


0.94 ±0.12 




060428B 


1 


1.70 ±0.14 


1.80 ±0.11 


2.35 ±0.25 


0.86 ±0.07 




060502A 


1 


0.45 ±0.05 


2.37 ±0.34 


1.21 ±0.19 


1.07 ±0.27 


1.510 


060502B 


3 


0.26 ±0.25 


0.95 ±0.55 




0.62 ±0.56 




060510A 


1 


0.76 ±0.08 


2.73 ±0.30 


1.05 ±0.07 


0.97 ±0.08 




060510B 


1 


0.78 ±0.05 


1.27 ±0.33 


1.57 ±0.35 


1.67 ±0.36 


4.900 


060512 


1 


1.43 ±0.22 


0.88 ±0.16 






0.443 
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Table 3 — Continued 



GRB 


D 


Pp 


(3 P d 


Pa 


Pad 


z 


060522 


1 


0.55 ±0.11 


0.95 ±0.13 


0.82 ±0.22 




5.110 


060526 


1 


1.00 ±0.15 


0.87 ±0.10 


0.89 ±0.14 


0.13 ±0.62 


3.210 


060604 


1 


0.88 ±0.30 


2.72 ±0.19 


1.36 ±0.20 


1.25 ±0.16 


2.680 


060605 


1 


0.30 ±0.11 


0.25 ±0.43 


0.89 ±0.09 


0.96 ±0.10 


3.800 


060607A 


2 


0.47 ±0.05 


0.94 ±0.10 


0.80 ±0.06 


0.86 ±0.09 


3.082 


060614 


1 


0.96 ±0.03 


0.88 ±0.12 


0.92 ±0.08 


0.77 ±0.12 


0.125 


060707 


1 


0.75 ±0.10 


0.84 ±0.24 


1.03 ±0.13 


1.20 ±0.46 


3.430 


060708 


1 


0.66 ±0.08 


2.22 ±0.51 


1.73 ±0.20 


1.47 ±0.16 




060712 


1 


0.82 ±0.20 


2.07 ±0.33 


1.36 ±0.34 


1.75 ±0.23 




060714 


1 


0.95 ±0.07 


1.36 ±0.31 


1.12 ±0.12 


1.36 ±0.27 


2.710 


060717 


1 


0.68 ±0.26 


1.23 ±0.42 








060719 


1 


1.03 ±0.08 


1.47 ±0.62 


1.36 ±0.11 


1.37 ±0.15 




060729 


1 


0.83 ±0.10 


1.64 ±0.05 


1.42 ±0.04 


1.39 ±0.05 


0.540 


060801 


3 


-0.55 ±0.18 


0.51 ±0.32 
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Table 4. Potential jet breaks. Jet break times, Tj, in seconds. For the top 8 GRBs the 
break must have occurred at or before the end of the plateau and the jet break time given 
is T a . For the bottom 4 the break was observed in the decay after the plateau and is listed 
in Table [2j ajd and (3jd are the indices measured after the break. The estimated jet angle 
9j (degrees) and jet energy log 10 (£' 7 ergs) are listed for GRBs with measured redshifts, z. 

GRB log 10 (7)) a jd Pjd z 9j log 10 (£ 7 ) 

050603 4.83 1^ 1.74 \f 2 0.70 ±0.10 

050730 4.13 tor 2.74 |;§f 0.62 ± 0.05 

050802 3.96 1.66 \f 4 0.81 ±0.09 

060413 4.58 £44 1.96 0.53 ±0.12 

060421 3.04 HI 1.25 \f 3 0.33 ± 0.37 

060526 3.84 |f 1.06 lf 9 0.13 ±0.62 3.210 1.5 49.4 

060605 4.16 1 ^ 2.04 ? : 39 0.96 ±0.10 3.800 2.2 49.3 

060614 5.00 1.98?;^ 0.77 ±0.12 0.125 10.8 49.5 



2.821 


3.0 


50.8 


3.970 


1.6 


49.9 


1.710 


2.2 


49.4 



050315 


5.48 


2-00 


1.29 ±0.31 


1.949 


6.8 


50.9 


050814 


4-931:1? 


1.80 3-16 


0.71 ±0.10 


5.300 


2.7 


50.7 


060105 


5 04 5 - 26 
d - u ^ 4.66 


3.56 HI 


1.07 ±0.03 








060607A 


A 81 63.51 
4 - 81 4.06 




0.86 ±0.09 


3.082 


3.4 


50.3 



Table 5. Comparison of the final X-ray decay index, a a d, with the optical decay index, 
a opt , for afterglows with no jet break. *This error is dominated by systematic curvature of 

the decay rather than photometric accuracy. 



GRB 


T a sees 


Otad (90%) 




050525A 


8.3 x 10 2 


1.30-1.48 


1.46 ±0.15 


050820A 


9.1 x 10 3 


1.14-1.21 


1.03 ±0.04 


060206 


7.2 x 10 3 


1.18-1.29 


1.26 ±0.2* 


060707 


3.8 x 10 3 


0.75-0.91 


0.97 ±0.09 


060729 


1.2 x 10 5 


1.19-1.34 


1.27 ±0.03 



